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We study the statistics of multi-photon absorption/emission processes in a mesoscopic ring 
threaded by an harmonic time-dependent flux For this sake, we demonstrate a useful analogy 

between the Keldysh quantum kinetic equation for the electrons distribution function and a Contin- 
uous Time Random Walk in energy space with corrections due to interference effects. Studying the 
probability to absorb/emit n quanta hu per scattering event, we explore the crossover between ultra- 
quantum/low-intensity limit and quasi-classical/high-intensity regime, and the role of multiphoton 
processes in driving it. 
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I. INTRODUCTION 

Recently a surge of interest in the dynamical proper- 
ties of mesoscopic/nanoscale electronic systems has mo- 
tivated a number of theoretical and experimental stud- 
ies on the physics of electronic devices subject to the 
driving of external fields. This main theme, pioneered 
in Ref. embraces a number of interesting issues such 
as the study of the influence of microwave driving on 
transport through chaotic scatterers^, the phenomenon 
of adiabatic quantum pumping- , as well as diffusion and 
localization^i^ii in energy space in quantum chaotic sys- 
tems/disordered quantum dots. In a broader context, 
the effect of the driving of external microwave fields has 
been shown to lead to an intriguing zero resistance state 
in quantum Hall systems-, and is currently studied as a 
tool to control the coherent dynamics of superconducting 
Josephson qubits^. 

Investigations of periodically driven mesoscopic sys- 
tems/quantum dots addressed mostly the limit of low 
intensity driving^i^ii. In this case, electrons have enough 
time to explore ergodically all available phase space be- 
fore performing a single photon assisted transition in en- 
ergy space. This makes it possible to use an effective 
time dependent Random Matrix Theory to describe the 
dynamics of the system. On the other hand, as beau- 
tifully shown by recent experiments in superconducting 
qubits^ and in the quantum Hall regime?, as the intensity 
of driving increases one should expect both an enhance- 
ment of the probability of single-photon processes, and 
the emergence of multi-photon processes/resonances in 
the dynamical properties of the system under study. 

The goal of this paper is to characterize the infiuence 
of multiphoton processes on the dynamical properties 
of mesoscopic electronic systems, concentrating on their 
effect on the electron dynamics in energy space (diffu- 
sion/localization). Diffusion and localization in energy 
space, as well as of multi-photon processes, have been the 
subject of a number of studies in the context of the optics 




FIG. 1: The physical system under study, a diffusive quasi one 
dimensional ring thread by a time dependent flux ${t). For an 
harmonic time dependence = $cos(a;i), the scattering 
of electrons off impurities induces transitions in energy space 
quantized in units huj. The statistics of such transitions and 
its physical consequences as the intensity of driving grows are 
described by Eq.® and Eq.lO. 



of complex atoms/moleculesiS. In these systems the un- 
derlying electron dynamics is typically very complex and 
a statistical description, either equivalent to random ma- 
trix theory or explicitly using it, is compulsory. In con- 
trast, in the present study we go beyond random matrix 
theory focusing on a model mesoscopic system, a diffusive 
quasi-one dimensional ring threaded by an oscillating flux 
(see Fig [Hand Eg. p^ ). were the underlying microscopic 
dynamics can be studied in detail. We explore how multi- 
photon processes proliferate as the driving amplitude is 
increased or the driving frequency is decreased, study 
the resulting crossover from ultra-quantum/low- intensity 
limit to quasi-classical/high-intensity limit, and extract 
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its physical consequences. On the theoretical side, we 
demonstrate and use extensively an interesting analogy 
between the quantum kinetic equation for the electron 
distribution function and the recursion relation defining 
a Continuous Time Random Walk^^ in energy space. 

The rest of the paper is organized as follows. In Sec.HIl 
we present qualitatively the results of our analysis of 
diffusion in energy space and of multiphoton processes 
based on the mapping of the problem onto a continuous 
time random walk in energy space. This mapping is de- 
rived in full detail in Sec. lIIII using the Kcldysh technique. 
Finally, in Sec. IIVI we present our conclusions. 



may immediately write 



t') ^'ilj{t-t'){S (n - / dt"evFE{t") ■h]) (3) 



where (*) = j dn (*) denotes the averaging over mo- 
mentum directions. For an harmonic time dependence 
E{t) = Eq cos{uit), and at low frequencies lot <C 1, one 
obtains 



Pn = / {Vn{t,t'))Tdt' = \20osin^ - ' - 



E, ( li^^ -1 



II. QUALITATIVE ANALYSIS 



In this section we start by summarizing the qualita- 
tive picture emerging from our analysis. The elementary 
time scale controlling the dynamics of energy absorb- 
tion/emission is the mean free time r. Indeed, in a dif- 
fusive quasi-ld ring threaded by a flux <i>(t) = $ cos(wt), 
energy changes quantized in units fko occur provided an 
electron scatters off an impurity. This is due to the fact 
that during the ballistic trajectory in between scatter- 
ing events the flux perturbation V{t) ~ —-^{1) v, v being 
the velocity operator, commutes with the unperturbed 
Hamiltonian Hq = therefore causing no transi- 

tions whatsoever. 

In the ultra-quantum limit of weak perturbations 
single-photon processes dominate. In other words, in one 
scattering event an electron may either absorb/emit one 
quantum, or scatter elastically. In particular, the prob- 
ability Pn to make a transition of energy O in energy 
space in one scattering event is given by 

Pn = (1 - p)S{n) + I [6iQ - hiu) + S{Q + hio)] , (1) 

where p cx ($/"I>o)^ <C 1, 'I'o being the flux quantum. 

On the other hand, in the opposite limit of high inten- 
sities of the perturbation it is natural to expect quasi- 
classical continuous energy absorption described by a 
Drude-like picture. According to this picture an elec- 
tron moving ballistically between two scattering events 
(at times t' and t, respectively) acquires an energy 

dt" evpEit") ■ h, where E{t) = -dtA{t) is the elec- 
tric field generated by the time dependent flux, and h is 
the momentum direction in the d-dimensional space. Let 
us introduce the probability density Vn{t, t') of changing 
the energy by between two successive scattering events 
at t' and t. Given the Poisson distribution 



tp{t-t') ^-e 

T 



\t-t'\/T 



(2) 



of time intervals \t — t'\, neglecting acceleration by an 
electric field E{t), and assuming isotropic scattering, one 



where (*)t = J^ji, ^(*) denotes averaging over the pe- 
riod, Ei{z) = /^°° e~^* ^ is the exponential integral func- 
tion, Kq{z) is the Bessel function, and = cEqVft. 

Multiphoton processes drive the crossover between 
discrete energy absorbtion in the ultra-quantum limit 
[Eq.([T])], and continuous energy absorbtion in the quasi- 
classical limit [Eq.Q]. As shown below, the crossover 
probability function Pq may be written as Pn — 
X^n Pn5{^ — hnuj). In the low frequency limit lot <C 1, 
and for isotropic 3d scattering, we obtain 



(5) 



where E — eEf)VpT / {hw), is a generalized hyper- 
geometric function, A,^ ^ (22'T[n + l/2])/(^(l -|- 
2n)r[n + 1]), a„ = {n + 1/2, n + 1/2, n -|- 1/2}, and 
&n = {n + 3/2,l + 2n}. These functions are plotted for 
selected values of n in Fig.l. 

At low intensities <C 1) the probability to ab- 
sorb/emit n photons is 



(6) 



Therefore, multi-photon processes are exponentially sup- 
pressed. Neglecting them we obtain Eq.(IT]) with p = 
£^/6. As the intensity grows, higher order processes be- 
come increasingly probable at the expense of single (or 
in general low) order ones, as indicated by the fact that 
for f > 2, Pi(£) starts decreasing. At f ^ 1 Eq.® can 
be approximated as follows: 



f \exp[— Ti/f], E <^n 



(7) 



Note that in the interval \ <SC_ n ^ E the probability 
of absorbing/emitting n photons decrease very slowly 
with increasing n which leads to a proliferation of multi- 
photon processes at large E . 

The two distinct regimes of rare (£ ^ 1) and pro- 
liferating (£ ^ 1) multiple photon processes have been 
discussed by Keldyshi^ in his seminal work on atom ion- 
ization. In this case the two regimes are classified by the 
ratio = i^t/w of the inverse time Wt = cEi^j \fml^ 
of tunneling through the potential barrier titled by the 
electric field Eq, Iq hui being the ionization threshold. 
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FIG. 2: The probability to absorb/emit n photons in a scat- 
tering event P„ plotted as a function of the intensity parame- 
ter S for selected values of n. At low intensities, single photon 
processes dominate. At higher intensities, higher order mul- 
tiphoton processes become increasingly important, driving a 
quantum-to-classical crossover [see text]. 



and of the frequency of the apphed electromagnetic field 
uj. The quahtative connection between Ref. 12] and our 
problem is obtained by identifying 7^^ with £, and /q 
with 1/t in the present analysis. 

Let us now make the qualitative considerations above 
more precise. As shown in Sec lIIII in terms of a Keldysh 
diagrammatic analysiSi^J^iH, the dynamics of energy ab- 
sorption/emission in the system at hand may be conve- 
niently described as random walk in the energy space de- 
scribed by the recursion relation for the electron energy 
distribution function ft{E) 



ME) 



dt' 



dnVn{t,t') ft'{E~n). (8) 



Neglecting weak localization effects controlled by the pa- 
rameter Xf/vpT <C 1 and effects of dynamic localization 
controlled by the parameter 6/ {vpTA{t)'^) <C 1 (where S 
is the mean separation of electron levels in a finite sys- 
tem)^, the kernel Vn, is given by the product of two 
functions Vn{t,t') = tp{t — t')pn{t,t'). The function 
tjj{t — t'), given by Eq.([21), is the distribution of the bal- 
listic time of flight \t — t'\, which may be interpreted 
as the continuous waiting time in between steps of a 
random walk in the energy space. The other function, 
po(i,t'), is the conditional probability to absorb/emit 
an energy during the ballistic flight and is given by 
pnit,t') = f dr]e-'"''prjit,t'), where 



(9) 



This result holds for a generic time dependence of A(i) = 
nx^{t)/L as long as e|A(t)|AF ^ 1. A random walk 
of the type defined by the recursion relation Eq.® is 
known in literature as a Continuous Time Random Walk 
(CTRW)". 

One may now easily derive the crossover probability 
function Eq.®. Indeed, in the case of harmonic flux 



Pn{t,t') 



5{VL — nhuj) Pn{t, t') 



(10) 



J2 



2 — fi ■ fix [cosiujt) — cosiujt')] 

UJT 



where J2n(,z) is a Bessel function and £ = eEoVpT /hw. 
Though pn{t,t') does not explicitly depend on r, the 
probability (averaged over the period T of flux oscilla- 
tions) 



Pr, 



dt'ij{t-t'){pn{t,t'))T, 



(11) 



of a multi-photon process between two successive scatter- 
ing events does depend on the mean free time r through 
the function ip{t). For wr ^ 1 one expands the dif- 
ference of cosines in Eq. ljlOp . introduce a new variable 
{t — t')/T and immediately concludes that P„ = Pn{£) 
is a function of £ — eEf)V ft / {fiio) . Performing explicit 
integrations in id leads finally to Eq.([5]). Similar results 
may be derived in the high-frequency limit ujt 3> 1. In 
this case where one may set ip{t—t') w 1 and observe that 
Pn = Pn{£) is a function of the r-independent parameter 

£ = cEqvf /^'^ ■ 

It is now possible to show directly that the discrete 
probability distribution P„ interpolates between ultra- 
quantum limit and quasi- classical continuous energy ab- 
sorbtion. First of all, one may directly compare Eqs.Q 
and (|4]) : calculating the average over one period in Eq. ([4|) 
(which is dominated by small t <^ uj^^ for f2 <C f^nd 
by ujt « 7r/2 for SI 3> rio) and replacing U, — > nhw one in- 
deed obtain the classical result of Eq. ([7]) • An alternative 
way to see the crossover, is to compute the moments of 
the number of absorbed/emitted photons ((?t-)™) exactly 
using Eq.®. One obtains 



(12) 
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The ultra-quantum limit corresponds to all ((n)^'") ~ 
i.e. keeping only the first term in Eq. (fT2|) . On the 
other hand, the classical distribution [3o?-case in Eq.(j3])] 
leads to moments ((Sl)^™) coinciding with the last terms 
in Eq. (jl2p upon the replacement n il/hw. This shows 
again that multi-photon processes drive the system at 
large intensities £ 3> 1 towards its quasi-classical limit. 

The table of moments Eq. ([12]) can in principle be ex- 
tracted from the smooth envelope of the distribution 
function f^^'"{E), using the standard techniques of the 
theory of random walksii to translate the properties of 
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FIG. 3: The self consistent Born approximation for the re- 
tarded/advanced Green's functions Ea. (|17p in its standard 
diagrammatic representation^ii^. The single/double hnes rep- 
resents the bare/full Green's function, wavy hnes correspond 
to the external drive V{t), while the dashed line is the disorder 
averaging. 



probability kernel into a complete characterization of 
its dynamics. It is clear that the finiteness of the sec- 
ond moment {(n)'^) implies to zeroth order a standard 
diffusive dynamics fr^iE) ^ EYfc[E/y/2D^ 

with diffusion constant in energy space De = " ^^^"^ ^ = 
Z3[ei?o]^, where D — Vpr/d is the diffusion constant in 
the real space. Higher moments, such as ((n)^), which 
in contrast to the second moment do contain information 
about the ultra-quantum to quasi-classical crossover, in- 
fluence higher order corrections. For example, up to first 
order in r/t we obtain 



1, 



-Erfc 



V2. 



A4 T e "2 
t 24V2^ 



z(z2 - 3) + . .(.13) 



where z = u/^/DEt and A4 = {{ny) / {{nYy — 3 is 
the kurtosis associated to distribution P„. The envelope 
ff™{u>) can in principle be measured by a tunnelling ex- 
perimen1)i5 . Note that although the term cx A4 in Eq. (fT3|) 
is a correction, it goes beyond the universal limit of RMT, 
corresponding to r — > 0, and is determined by the details 
of the semiclassical electron dynamics (e.g., smooth dis- 
order/anisotropic scattering, quantum isotropic scatter- 
ing). 



III. FORMAL DERIVATION 



Let us now outline the formal derivation of the map- 
ping of the dynamics of the distribution function onto a 
continuous time random walk in energy space, Eq.®. 
We adopt a model of free electrons in a Gaussian 6- 
correlated static impurity potential U (r) (which corre- 
sponds to isotropic scattering amplitude) coupled to an 
external time-dependent vector potential A{t), through 
V{t) = ~evA{t). The Hamiltonian takes the form 



Zm 



V{t). (14) 
- r'), v being the den- 



where {U{r)U{r')) = l/{2Tn'T)d{r- 
sity of states at the Fermi level. 

The distribution function ft{E) may be expressed as 
ft{E) = 5 [1 — / drje'"^^^ ht{ri)] , where in turn /it (77) can 
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FIG. 4: (a) The contribution to the Keldysh Green's func- 
tion 5G^. Eq. (|19p in its diagrammatic representation to low- 
est order in the small parameter 5/{vpTA{t)^). The zig-zag 
line represents the insertion of ho{ti — t2)[V{t2) — V{t-i_)]. (b) 
The diagrammatic representation of the diffuson appearing in 
5G^. (c) The perturbative expansion of the triangular vertex 
in the external drive V{t). 



be written in terms of the disorder averaged Keldysh 
Green's function^'^^ as 



htiv) 



2ttv 



dr 

v 



{G^{t + r,/2,t~r^/2;v,v)).il5) 



Let us now exploit the structure of the perturbative 
expansion of in the time dependent perturbation 
V{t). As shown in detail in RefQ, the noninteract- 
ing nature of this problem makes it possible to iden- 
tify two contributions to the Keldish Green's functions, 
(G^) =G{f + SG^. The first 



G, 



K 



dt"{G^it,t',r, r))hoit" -t') 
ho{t-t"){G^{t",t',r,r)), (16) 



represents physically the unperturbed distribution func- 
tion. Indeed, calculating the disorder averaged retarded 
and advanced Green's functions (G'''°), within the self 
consistent Born approximation [see Fig. ([3])] , one obtains 

(G'''''(i,t'))p = TiOi±tTt')e-"'''^'~*">e^^ 

^-ivf-n-f*, dt"eA{t") ^ j-j^y-j 

which immediately implies {G^''^{t,t' ,r,r)) = ^iTri'S{t — 
t'). Therefore, G^(t -h 77/2, i - r//2, r, r) ^ -2nmhoi'n). 
The second contribution to G^, 

SG^ = J dtidt2dri(G^{t,ti,r,n)ho{ti-t2)[V{t2) 

-V{t,)]G-{t2,t\rur)), (18) 

describes energy absorbtion from the time dependent 
field. Performing now the disorder average in Eq. (jl8p 



5 



the product of a retarded and advanced Green's function 
appearing in Eg. p^ generates a diffusion propagator. 
More specifically, SG^ admits a diagrammatic represen- 
tation in terms of a loose diffuson [see Fig.(|3]-(a))], which 
formally amounts to the equation 

SG'^ = 2Triiy [ dt' dt"V^{t,t')C^{t' ,t") ho{r]), (19) 



where, D is the standard diffuson [see Fig.(j4]-(b))] solu- 
tion of the equation 

D-^®Dr, = Dr,{t,t')~ j dt"U^{t,t')D,j{t",t') 

= S{t-t'). (20) 

Neglecting interference effects, the kernel 11, as well as 
the vertex C are given by 

U,it,t') = J drfTr{G^iU,t'+))(G-it'_,t^))/{27T,yT) 

C^{t,t') - J d7j'Tr{G%U,t'+)) {G-{t'_,t^))[V{t'_)) 

where i± ^t±ri/2 and 4 ^t'±r]'/2, [see Fig.(|4i;b)-(c))]. 
The Tr symbol stands for the trace over the coordinate 
indices; in particular it implies J dp — J de{p) J dh in the 
momentum representation where the disorder averaged 
Green's functions (G"''°(i, t'; p)) are diagonal. Preform- 
ing explicitly the trace, one obtains 

Ur,{t,t') = e{t-t')^b{t~t')pr,{t,t'), (21) 

Cri{t,t') = e{t-t')ij{t--t')dt'PrM.t'). (22) 



where ip{t) is given by Eq.((2|) and 
Pri{t,t') — J dn exp 



t+ri/2 

ivpn-Uxl I dt"eA{t") 

t'+r,/2 



t-ri/2 

dt"eA{t") 

't'-r,/2 

Finally we may express ht{ri) in terms of 2?, C as 



(23) 



htiv) ^ (^l- J dt'dt"V„{t,t')C^{t',t")^ hoiv). (24) 

Let us now show that the distribution function is in 
the kernel of the inverse diffusion propagator, i.e. 



D~^®ht{Ti) = 0. 
First of all notice that 



(25) 



= -[V-']^{t,t')+C^{t,t'). (26) 



n 



R A/\ R 

+ 



FIG. 5: Loop expansion for Pn{t,t'), the insertion C being a 
cooperon^. 



Now acting on Eq. ([24|) with the operator D ^ we obtain 



dt'Crjit,t')-T / dt'dt'Il,j{t,t') 



dt'£rj{t,t') 



^ 0. 



(27) 



Notice at this point that Eq. ([25|l is equivalent to the re- 
cursion relation 



(28) 



bmce J dt'Il^^o{t,t') = 1, the latter may be equivalently 
stated as 



dt' n^(t,i') fr{vi) 
t 

dt'V„{t,t')ft,{i^), 



(29) 



where Pjj{t, t') = Tp{t—t')pr^{t, t'). The Fourier Transform 
of this equation with respect to rj gives the recursion re- 
lation defining the continuous time random walk, Eq.®. 

It is natural at this point to ask whether the formal de- 
scription of the time evolution of the distribution func- 
tion as an effective random walk may include interfer- 
ence/localization effects as well. Indeed, performing a 
one-loop analysis with accuracy £^ one may shown that 
the structure of Eq. ^ persists with a probability kernel 
schematically represented by the diagrams in Fig. ([5]) . It 
is clear however that, in the presence of interference, the 
Markovian nature of the random walk is lostii. In partic- 
ular, upon time integration we obtain a probability dis- 
tribution Pji given by Eq.© with p{t) ~ ^ {l - , 
where the driving time t is the time since the turning 
on of the perturbation, t* — 2tt^ D e / {^"^ S"^) is the lo- 
calization time in energy space^, and we neglected all 
corrections independent of t. This result amounts to 
the weak localization suppression of the absorbtion rate 
W{t)/WQ ~ 1 — \/t/t* due to weak dynamical localiza- 
tion^. 
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IV. CONCLUSIONS 



In conclusion, we have studied the problem of en- 
ergy absorbtion/emission in a mesoscopic ring threaded 
by an oscillating flux, focusing on the influence of mul- 
tiphoton processes, and on the multiphoton driven a 
crossover from ultra-quantum/low-intensity limit, and 
quasi-classical/high intensity regime. We have shown 
that the dynamics of the distribution function may be 
mapped onto a continuous time random walk in energy 
space. Though in the present paper we focused on the 



effect of a classical driving, recent advances in the field 
of circuit QED^ strongly suggest the possibility to in- 
vestigate the role of single and multi-photon processes in 
the case of quantum driving, an interesting problem that 
remains a challenge for future work. 
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